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Proposed for use in problems of orientation, navigation and Euler rigid body dynamics (RB) conjugate ("tangent" 
and "cotangent") vectors of finite Euler rotation. Modules of these vectors are proportional to tangent and cotangent of a 
quarter angle of rotation. Contrary to the classical Rodrigues ' vector, they do not become infinite within the complete RB 
rotation. Kinematic differential equations for these vectors have common solution (in the Cauchy's form). New dynamic 
polar differential equations are obtained from the classical dynamic Euler-Poisson equations as a result of replacement of 
angular velocity vector by the vector angular momentum, as well as replacement of three directional Poisson cosines by the 
conjugate vectors coordinates. In the case of Euler, three obtained scalar dynamic differential equations have two classical 
first integrals (of energy and areas). Three new dynamic equations are considered as integrable ones. These equations 
explicitly determine the whole RB orientation, not just its vertical, like three Poisson equations do. Conjugate vectors X , p 
are considering as parts of vectors with non-traditional non-normalized (non-Hamiltonian) rotation quaternion and 
biquaternion. Such non-normalized quaternions can be effectively used instead of classical normalized quaternions and 
biquaternions (with parameters of Rodrigues-Hamilton) in different «particular» tasks ofinertial orientation and navigation, 
including RB orientation control problem. 

1 . In inertial orientation and navigation are of particular interest dimensional rotation vectors having a 
minimum number of independent parameters, corresponding to three degrees of freedom of the RB with a fixed 
point [1]. In classical Euler' s problem of RB dynamics with a fixed point [1], the differential Euler-Poisson 
equations are commonly used. Three Poisson equations determine only the earth vertical reference in relation to 
RB, not the whole RB orientation. Moreover, it is presumed possible to determine RB orientation through the 
Euler angles with additional integration of expression for time derivative of RB precession angle [1, page 16]. But 
the direct integration of only three Poisson equations basically does not solve the problem of determining RB 
orientation in three-dimensional space. In order to explicitly determine RB orientation, it is necessary to have at 
least six directional cosines. However, in this case, one must integrate nine differential equations instead of six 
Euler-Poisson equations. In this regard, the possibility of using well-known [2-10] three-dimensional rotation 
vectors in RB dynamics is appealing, because they have three independent coordinates. 

2. Classical Rodrigues's [2,3] vector has a module with tangent tan(cp/2) of half-angle (p of finite 
rotation of RB. This vector can not be used in the condition of (p = n in the problems of RB dimensional 

orientation. Conjugate rotation vectors x = xk , p = pk , where x = k x tan(cp / 4) , p = k p cot(cp / 4) ( k z , k p are 

arbitrary constant coefficients), k - unit vector of Euler axis of RB finite rotation are of the most interest (by 
0 < cp < 2n ) because of its conjugate features [3]. 

The article deals with the conjugative polar non-linear vector kinematic differential equations of RB 
rotation in the transformation of lineal vector a>(?) of RB angular velocity [3,7]: 

f* = t'B (t) &, p = p'Q T lp) a, (l) 

where x - dx I dt , p = dpi dt - local derivative of vector in the time t (derivative towards some 
connected with RB coordinate base); x'=dx/dq>, p'=dp/dq> - partial derivatives of modules i, p of vector 

of angle cp, defined as function x' = (k x + x 2 1 k T )/4 , p'= -(k p +p 2 /£ p )/4 ; 9 (r) , 6^ p) - orthogonal operators: 
e (r) =E + 2(fc r T + T 2 )/fe + r 2 ), 
Q T {p) =E + 2(-k p R + R 2 )/(k 2 p+ p 2 ). 

Therein T , R are skew-symmetric operators of vector multiplication [3,5,8], satisfies the identities 
Tx = xxx = 0, Tp = pxp = 0, 0 - zero vector, E - unit operator. Operators 0 (r) , 6 (p) satisfy the identity 

(0 (r) 0 (p) ) = E , determined conjugacy [3] of vectors x , p and conjugate equations as features of "duality" and 
isomorphic correspondence. 



3. Conjugate equations have the first («trigonometric») integral in the form of scalar product 
(x • p) = xp = k x k p = C (arbitrary constant) and admit simple and graphic kinematic interpretation - polar 

wobbling-nutational model [3,9] with arbitrary vector ro(f) and for any RB. In the first equation vector a> is 

transformed (at every moment of the travelling time of RB) by operator 0^ into vector ro T as a result of 

precession b ro on the angle \|/g T =q>/2 (rotation vector S with module ra(f) about the surface of some circular 

cone of «velocity» precession cone angle 2u ro ) around Euler axis with unit vector k . Angle v m - nutation angle 

(vector deviation ro from the unit vector k ). And then a vector a> x is multiplied by scalar operator z'E . At the 

second equation vector © precesses (rotates n the surface of that) at the angle v|/ ep = (n + q> / 2) and is multiplied 

by scalar operator p'E . Angle u ro of nutation is determined from scalar product (ro-x) , (x*-x) or (ro-p) , 

(P -P). 

On the basis of such interpretation - equation (1) models come out, for example, kinematic 
nonholonomic equality: 

((5 -b) 2 +{{dxb)-{dx b)))/((5 ■ a%b ■ b)) = sin 2 ^ + cos 2 ^ = 1 , (2) 
where a = (ox x , b = (0 (r) <»)x f . 

The equality with arbitrary constant C T =k T , comes out from equality, for example: 
2 (S-a) + (a-b) 

( -T =T !T~^~ ' (A) 

(a -a) -(a-b) 

Equalities, analogous to that equalities, comes out of vector change x to vector p . 

4. Under the numerical integration [2-4] of conjugate equations such values are to exist k x * and k p * of 



coefficients k x , k p , where equality is provided x'=l h p'=l . Then, for example, k x *=2 + ^<\-% 2 , 

4-p ,(x<2,p<2) vector transformation ro into derivatives x , p will be exercised by 

orthogonal operators 0 r » , Q T p * , where at every step of calculation new values of "adaptive" coefficients k T * , k p * 
(instead of constant coefficients k x , k p ) are used. But in fact integrate conjugate equations in the form of orthogonal 

transformation T* = V T ,G), p* = V p *G) when known (approximately «measured» by integrating gyroscopes) of 

kinematic integrals (vector modules x ,p ) x*«p*«g + c 9 «(p from scalar functions dx* I dt = x*(f) , p*(f) , 
ro(f) (c q - arbitrary constant, ro(f) - module of vector). 

5. Conjugate equations have general solutions in the Cauchy form [3,6]. General solution of the first 
equation in (1) determines the form of group operations of multiplication and division in group of Lie vectors x 

[3,7]. In this group inverse element x _1 is equal to adverse vector, i. e. x _1 = -x , but unit element is equal to zero 

vector 0 . 

6. New polar dynamical differential equations of RB rotation are obtained on the basis of kinematic 
conjugate equations. 

In the classical Euler case, for example, the problem of integration of the system of six dynamical 
differential equations of Euler-Poisson [1] reduces [7,9] to integration of the system of three dynamical equations 
with two independent classical first integrals (the energy and area integrals). These two integrals are enough to 
consider the system of three equations to be integrated [1]. 

In this case change is introduced in the conjugate equations ra = S - g , where g - constant vector of 

angular momentum (constant module and in the line of supporting basis I ); S~' - inverse operator S of RB 
inertia (constant in connected with RB basis J ). 

Contained of conjugate equations, for example, polar matrix differential equation with vector coordinates 
x has the form [9]: 

^J=A Q (T)jSj%)j)Q T (T )j8i, (4) 



where Tj = [x ; i Tj 2 ^nf - column matrix with derivative coordinates x in basis J ; Qmj - 
matrix ( 3 x 3 ) of operator 0^ r) in basis J ; Sj - diagonal matrix (3x3) with three constant the main moments of 

RB inertia; S] 1 - inverse matrix, gj = [g n g j2 g^Y ~ constant column matrix with vector coordinates g in 
basis I . 

Three vector coordinates x definitely determine RB orientation (in contrast to three direction cosines 
determined orientation of just one unit vector of supporting basis I in Poisson equations). In the problem of 
synthesis of RB orientation control law (determination of control moments) applying functions of Lyapunov 
equation is used together with classical Euler's dynamical equations. 

7. Conjugate vectors x , p are considering as parts of vectors with non- traditional [10, 11] non- 
normalized (non-Hamiltonian) rotation quaternion and biquaternion. Such non-normalized quaternions can be 
effectively used instead of classical normalized quaternions and biquaternions [2, 10] (with parameters of 
Rodrigues-Hamilton) in different «particular» tasks of inertial orientation and navigation, including RB 
orientation control problem. 
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